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Abstract 

■ Equations for the center of gravity of the shower originated by high energy proton in the atmosphere 

are written and, within certain simplifications, solved for the case of logarithmically decreasing interaction 
length of hadrons in the air. Obtained expression provides transparent view of the way in which hadronic 
interaction characteristics determine the longitudinal shower development. 

8 ■ 1 Introduction 

Since long ago, numerous attempts have been undertaken to explicitely connect air shower longitudinal profile, 
in particular the shower maximum depth, with hadronic interaction characteristics. 

Very approximate approaches were tried, from toy models to extensions of Hcitlcr model for the electro- 
J> | magnetic shower Q] to the hadronic shower [2j [3l [4] and all these have proven to be of not big quantitative 
help. 

A direct way to establish the connection is use of the cascade theory. A problem is that the shower maximum 
is an inconvenient quantity for treatment by cascade equations. Rather, a convenient quantity is the shower 
center of gravity (CG). It differs from the shower maximum basically by a shift. Thus, with CG one can trace 
the elongation rate of the shower maximum. Even better suited CG for studying the difference in the elongation 
^ rate between shower simulation codes using different interaction models and for separating main factors in the 

interaction properties determining this difference. 

In this work, equations are written for the center of gravity of the proton initiated shower and within certain 
simplifications an analytical solution is derived. 

2 Derivation 

We consider dependence on energy of the shower center of gravity: 
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X(E) = J X N(X) dX/ J N(X) dX = J X N(X) dX/ 
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For convenience, further we measure depth in the radiation units, i.e pass to variable t — X/ Xq. Also, we 
measure energy in units of the critical energy E c . 

The proton primary is considered. In simplifying assumptions of Feynman scaling and neglecting production 
of other than nucleons and pions particles and neglecting pion decay, the system of equations for the nominator 

oo 

$ = J t N(t) dt looks: 
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*o(-E) = 2 / — <f>^(Ex)dx 
J x 



Here N denotes nucleon, 7r denotes charged pion, denotes 7r°, 7 denotes photon and dn^j/dx defines dis- 
tribution over Feynman x of the secondary particle of type j produced by the primary particle of type i. 

For the center of gravity of the electromagnetic shower, ingnoring the Landau-Pomeranchuk effect, we can 
obtain: 

t-y(E) = In E + S 

Exact value of 5 is unimportant for our goals, in the approximation B of cascade theory it proves to be close 
to S = 1.7. 

Accordingly 

<$>-y(E) = E (In E + S) 

Let's apply Mellin transform: 



We obtain: 



F(s) = J E- S ~ 1 F(E) dE 
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Solution of the system of equations is: 



*jv(s) = 
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Here fij(s) — [ x s — ^—^-(x)dx 
Jo dx 



Let's assume logarithmic dependence of the interaction length on energy: 

\{E) = A {1 -alogCE)} 

In that case 

A(s + 1) = A ' 



s-1 (s-1) 2 

Let's apply inverse Mellin transformation to the solution: 

& N (E) - — <j> E s ® N {s)ds 
2-Tvi J 

<J?iv(s) is a sum of two terms with poles at s=l: 



s - 1 n v ' (s - l) 2 



*£>(•) = 



(a + l)(l-/jvjv(«)) 



(a + l)(l-/jvjv(«)) 



/jVtt(s) /tto(s) 

— — 7 \ h Jtvo(s) 

1 - f-xn(s) 



, Affl | /tv.(s) AW 



1 - fiviv(s) 



+ /jvo(s) 



1 - /jvjv(s) 1 - /jvjv(s) 
A^ } • ajv /jvtt(s) A^ • a,,. 



1 — /iVTv(s) 1 — /jVJv(s) 1 — /iT7r(s) 



The complex integral equals to sum of residues at these poles. 

— ^ / E s & N (s) ds — Sres $jv 
The contribution of the first pole is 

= i) = B (5 + ; f Ajv + — A,, , r . 

I 1 - 9nn V 9^0 / J 

where <7ij = = 1) and we have taken into account that .Qtvtv + Qjvtt + Sato = 1 and + g^o = 1. 

The contribution of the second pole is 

lim^{^)(3)}, 

which equals to 

E { InE — — + f' NN ll£lL f' N0 _|_ . ill® ^EZE 
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where 
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Let's denote 
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where index X in the inclusive distributions implies production of particles of any kind. 
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We obtain for the center of gravity 
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This is equivalent to the interaction lengths being taken at some effective, reduced relative to the primary 
ones, energies: 

E e J f = El exp ( f ' NN ) and E e ff = E e J f / exp ( ^JL + Ik.) 

\i-g N Nj 1 ' Vfl/VTr g-KoJ 

Finally, moving back to depth in g/ cm 2 and energy in GeV: 
X N (E) = X (in ^ + 8 - I) + ( X N (E e J f ) + X ■ » N + [\„(E e Jf) + X ■ ^] } 

\ E c 2 J 1 — (/TV TV I g-rvO J 



Interaction lengths are expressed trough inelastic cross-sections as A 



where A is atomic mass 



and Na is Avogadro number. 
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3 Discussion 



The final expression for the shower center of gravity explicitly splits into two terms: the center of gravity of 
the purely electromagnetic cascade at primary energy and a modification of this by the hadronic cascade. The 
latter is determined by two competing oppositely directed processes: i) carrying through energy by hadrons, that 
elongates the shower, ii) energy dissipation in the hadronic interactions due to which electromagnetic subshowers 
start at smaller energies, and because of the logarithmic energy dependency of their center of gravity that results 
in shortening of the total shower. The first proccess is represented by A terms in the final expression, the second 
one is represented by n terms. 

Three main quantity characterizing hadronic interactions are usually considered as governing longitudinal 
shower development: inelastic cross-sections, inelasticities and mean multiplicities. Cross-sections directly enter 
the final expression for CG through interaction lenghs A. Inelasticity, or more generally relative energy transfers 
between different particle types, also directly enter the expression as gij integrals (inelasticity is the denominator 
in front of the curly braces). Whereas multiplicity does not enter the expression as such. Instead, energy splitting 
is represented by the integral over inclusive distribution with additional weight xlnx relative to the integral for 
the multiplicity. The meaning of this weight is clear: each produced particle contributes to <fr with the weight 
which is the product of its energy and of the center of gravity, which is proportional to the logarithm of this 
energy. 

The obtained expression can provide only semi-quantitative results because of made simplifications, most 
important of which are neglecting Feynman scaling violation, production other than nucleons and pions particles 
and decay of charged pions. Nevertheless, difference in the elongation rate predicted by different generators 
should not be severely influenced by these simplifications. 

4 Conclusions 

Equations for the center of gravity of the shower from primary high energy proton in the atmosphere are written 
and, within assumptions of Feynman scaling, neglecting production of other than nucleons and pions particles 
and neglecting charged pion decay, solved for the case of logarithmically decreasing interaction length of nucleons 
and pions. Hadronic interactions are represented by inelastic cross-sections, taken at some effective energies, 
and integrals over inclusive distributions of two types: relative energy transfers between different particle types, 
like inelasticity, defining elongation of the shower via hadronic cascading, and integrals with specific weight 
xlnx, reflecting energy dissipation and leading to shortening of the total shower. 
Results of applying the obtained expression will be presented in next publications. 
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